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Abstract 



We measure field-strength and their correlators in presence of a static qq pair by 
numerical simulations. We give an interpretation of these data in terms of quadratic 
and quartic cumulants. 



1 Introduction 



A study of the field-strength distribution between static quark and antiquark yields 
a detailed picture of the string formation and helps to understand the mechanism of 
confinement. In particular, comparing the field distributions with those of the dual 
Meissner model one can clarify the viability of the popular confinement mechanism. 

In non-abelian theory there are two ways of measuring the field-strength distribu- 
tion: using connected p c or disconnected p dlsc plaquette averages around the Wilson 
loop [1]. While both reduce to the same quantity in the Abelian case, in the non- 
abelian case two measurements yield independent information. Recently [1] p c and 
pdisc h ave been measured for some components of the field F^ using the cooling 
method [2]. While the signal for p dlsc was too small, the distribution of E n off the 
string axis ("the string profile") and along the string from p c was found with good 
accuracy. Another important quantity measured in [1] was a double plaquette corre- 
lator 7 c (x, x'), yielding additional information on correlation of field strength at two 
different points x, x' of the string. 

The purpose of the present paper is: i) to extend the Monte-Carlo (MC) mea- 
surements to all orientations of the plaquettes in p c and 7 c (rc, x'), thus yielding the 
most complete information on the field distribution in the qq string; ii) to give inter- 
pretation of the results in terms of simpler quantities - field strength correlators like 
(F llu (x)F pa (y)) (cumulants). 

Using cluster expansion and non-abelian Stokes theorem [3], one can express p c 
and 7 C in terms of cumulants and keep the lowest ones (physical arguments in favour 
of dominance of lowest cumulants are given in the review paper [4]). The latter have 
been recently measured on the lattice [5], yielding a rather small correlation length 
T g = 0.2 fm. Now, using the cumulants, one unambiguously predicts p c and can 
compare it with MC measurements. This comparison with data from [1] was done 
in [6] and a good agreement of measured and calculated string profile was found. In 
particular it was clarified how a small correlation length T g = 0.2 fm yields a bigger 
radius of the string - 0.5 fm. 

After presenting new measurements in sections 2 and 3, we make in section 4 a 
more precise and extended calculation of p c in terms of the bilocal cumulant and 
compare it with MC data. In particular we predict vanishing of some p c (F fJ-u ) on 
symmetry grounds and observe it explicitly in data. We perform the same analysis 
for 7 C and predict vanishing of most combinations of plaquette orientations. The 
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dominating structure - j c (E\\,E\\) - is expressed in terms of the quartic cumulant, 
which has never been measured on the lattice; the data for 7 C allow to make some 
estimates of it. 

A short summary of results is given in conclusion. 

2 MC study of the field strength tensor. 

Using MC technique, we study the spatial distribution of the components of the field 
strength tensor in presence of a qq pair. This generalizes to all the components of the 
field the results already obtained in [1]. Following [1], we define: 

c {tr{WLP, u {x h x i _)L+)) 
P » v (TrW) { ' 

where W is a Wilson loop, L is a Schwinger line and is the part of the plaquette 
proportional to the a- matrices, oriented in order to give the desired component of 
the field. The coordinates measure resp. the distance from the edge of the 

Wilson loop and from the plane defined by the loop, as shown in Fig. 1. In the naive 
continuum limit a — > 0, 

P%, - a 2 {F, u ) qq (2) 

We have used a 16 4 lattice, taking a 8 x 8 Wilson loop and /3 = 2.50, which is inside 
the scaling window for the fields. Moving the plaquette in and outside the plane 
defined by the Wilson loop, we obtain a map of the spatial structure of the field as a 
function of x\\ and x±. Using a controlled cooling technique (see [1,2] and references 
therein), we eliminate the short-range fluctuations. The long-range non-perturbative 
effects survive longer to the cooling procedure, showing a plateau of 10-14 cooling 
steps, while the error becomes smaller. A similar behaviour has been observed for 
the string tension. The cooling technique allows us to disentangle the signal from 
the quantum noise with a relatively small statistics. The general patterns of the field 
configurations are briefly resumed in the following figures. Figure 2 represent a detailed 
map of the spatial behaviour of the longitudinal component of the chromoelectric field. 

• Varying x\\ at fixed x±, we investigate the structure of the fields in the direction 
of the axis joining the qq pair. In Fig. 3, we show E» and Bj_ as functions of x» 
for x± = 0, i.e. on the qq axis. We find that the electric field remains constant 
and magnetic transverse field vanishes, as expected on symmetry grounds. 
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• Varying x± at fixed x\\, we study the transverse shape of the fields. Fig. 4 
illustrates the behaviour of the E\\ component vs. x±, for different values of xf 
the field remains constant with respect to x\\ also outside the plane defined by 
the Wilson loop, as long as we remain inside the string (i.e. for x\\ = 3,5). A 
detailed study of the transverse shape will be given below. 

We find that the parallel electric field is squeezed in flux tubes, as already found 
in [1]. The results in [1] were consistent with a gaussian behaviour of the flux tube 
profile inside the string. In order to estimate the size of our tubes and to check the 
consistency of the result with previous measurements, we have again performed a fit 
of the transverse shape inside the string (for x\\ = 3) with the function 

E\\ = exp(« - n 2 x\) (3) 

finding \x = 0.30±0.01, with x 2 /d.o.f. = 0.993. This indicates that the flux tubes have 
a transverse size of the order of 3 lattice spacings at f3 — 2.50, which corresponds to a 
physical value n phys /Ai a tt = 85 ± 4, in agreement with [1]. In what follows, motivated 
by the measured form of the field strength correlators [5] and by the analysis in terms 
of cumulants [3], we will find that data are equally consistent with the form eq. (18), 
which is also suggested by the mechanism of confinement via dual superconductivity. 



3 Field strength correlators 

In the last years, a systematic study on non-perturbative effects in QCD in terms of 
the gluon field strength correlators has been developed (see Ref. [3] and references 
therein) and the behaviour of these correlators in the vacuum has been investigated 
by lattice simulations [5]. As pointed out by the authors of [3], studying the field 
correlators in presence of qq pair, could provide further informations for describing 
color confinement. Therefore we measure the operator: 

where 

V PP , = P» u LP> a L + - \p» v TrP' pa (5) 

where W is a Wilson loop, S is a Schwinger line connecting the Wilson loop to the 
Vppi operator, P and P' are two plaquettes, located at x and x' respectively, and L 
is a Schwinger line connecting them. 
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In the naive continuum limit, we have 



1 c = a 2 (F')+a 4 [(FF') qq -(FF') ] 



(6) 



where F and F' are respectively the field components at x and x'. 

Varying the orientations of the two plaquettes, we obtain the different components 
of the correlators. The measurements have been done on a 12 4 lattice at (3 = 2.50 
using 6x6 Wilson loop. Again, we used controlled cooling to reduce the fluctuations. 
We have measured 7 C with the following two types of orientations of F and F'. 

(i) In the first case: 

• P is held fixed on the qq axis at 1 lattice spacing from the border of the Wilson 
loop, while its orientation is varied in the 6 possible directions; 

• P' is moved in- and outside the plane of the Wilson loop, its orientation is kept 
fixed in the En direction; 

• x\\ and x± identify the position of P' with respect to P. 

(ii) In the second type of measurements: 

• both the position and the orientation of P are kept fixed; the plaquette is in the 
same position as before and its orientation corresponds to the E\\ component; 

• P' is moved as before and its orientation is changed. 
We finally define the irreducible correlator 7 C as follows 



From (7) it is clear, that 7 C contains only double plaquette correlations. 
Most of the data for 7 C and 7 C are compatible with zero net effect, within two standard 
deviations. In Table 1, we report the data for 7 C in case when both plaquettes P and 
P' are kept fixed in the En direction. In the next Section we compare the Monte-Carlo 
measurements with the predictions from the cumulant (cluster) expansion and will 
see that indeed all orientations except E\\ , E\\ should give zero result due to simple 
symmetry arguments. 

4 Extracting bilocal and quartic field-strength cor- 
relators from the Monte— Carlo data 

For the contour C shown in Fig. 1 we denote direction along the qq axis xn = x±, 
while that of x± — X2, and the Euclidean temporal axis is x±. All the construction in 
Fig. 1 is taken at a fixed value of x^. 



(7) 
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Using the non-abelian Stokes theorem and the cluster expansion [3] for p° in (1) 
one has (see [6] for details of derivation) 

p c fll/ (x 1 ,x 2 , x 4 ) = a 2 J dauix'^x'^A^ (8) 

where 

A^ = Ytr(E l {x[,0,x' 4 )<5>F IIU {x 1 ,x 2 ,x 4 )<5> + ) + (9) 

$ is the parallel transporter (Schwinger line) from the point (x[, 0, x' 4 ) to ( 

and dots imply contribution of higher order cumulants, containing additional powers 

of Ex [3]. 

We shall keep throughout this Section only the lowest cumulants (containing low- 
est power of Ei) and compare our prediction with the MC data of previous sections. 
The bilocal correlator A^ u can be expressed in terms of two independent Lorentz 
scalar functions -D((a^ — x'^) 2 ), Di^x^ — x' ) 2 ) (see [3] and the appendix 1 of the last 
ref. in [3]) 

A 14 = D + D 1 + {h\ + h 2 )^ (10) 

A 24 = (h 1 h 2 )^- , A34 = (Ms)^f (11) 

, , dDi . , , dDi . . 

A23 = 0, A13 = hzh 4 — ; A12 = ( 12 ) 

Here = (x - z%. 

Since all construction in Fig. 1 is at x 3 = x' 3 = we have h 3 = and hence 

P23 = P34 = Pn = (13) 

The only nonzero components are A 14 ,A 2 4 and A 12 . For the latter the contribution 
to p c can be written as 

Pi 2 (xi,x 2 , x 4 ) = a 2 dx\ J 2 t dx' 4 (+x 2 )(x 4 - x i) ~^ ' ( i4 ) 

When X4 = ( and this is where measurements of p\ 2 have been done), p\ 2 vanishes 
because of antisymmetry of the integrand in (14). 

Hence only p\ A and p 24 are nonzero, and only those have been measured to be nonzero. 

To make comparison with data more quantitative, let us exploit recent MC cal- 
culation of D and Di [5], which imply that both D and Di are of exponential form 

D^h 2 ) = £>i(0) exp(-f/i/i); D(h 2 ) = D(0) exp(-ph) (15) 
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ZMO) « -£>(0); fjL^fjL 
Inserting this into (8), (14) we have 

P c u(x 1 ,x 2 ;0) = a 2 J q R dx[ J T z dx' 4 [D(0) + ^(0) - M±Ml £ l(0 )] e -^ (i 6 ) 

with 

hi = —x' A , hi — Xi — x[ , h 2 = h\ + h\ + x\\ 
For similarly one obtains 

p c 2A (x u x 2 ;0) = -a 2 px 2 [* dx[ f~ dx 1 ^ 1 ~ X ' l) D^e^ (17) 

Jo J-\ 2h 

From (16) and (17) one can deduce, that 

(i) P24 should vanish for x 2 = 

(ii) P24 changes sign for X\ — ^, i.e. in the middle of the string length, 
(hi) P24 is about 1/3 of p^- 

All properties (i)-(iii) are supported by the data. 

Finally, we can make a detailed comparison of our prediction for p\ A in (16) with 
data. One obtains a simple analytic result for p\ A [x 2 = x±_) in case of a very long 
string. The transverse shape measured at the middle is given by [6] 

Pu = ^W)(l + px 2 ) - DM\i^ 2 ) 2 ]e~^ (18) 

As shown in [6], this shape is in good agreement with the previous data, obtained in 
[5]. Here we calculate p\ A as a function of x\,x 2 from (16) keeping -Di(O) = |D(0). 
We then fit the data for x\\ = 3 to evaluate p and a 2 D(0). We find: 

jf«0.19fm, a 2 D(0) w 3.92 x 10 7 

with a x 2 /d.o.f = 0.17. 

The value of p is in good agreement with [5], while we find that a 2 D(0) is one order 
of magnitude smaller than in the previous measurements. We recall that our data 
here are obtained for SU(2), while in [5] the gauge group was SU(3). This should 
account for the one order of magnitude between the two results. 
These results allow to predict all curves for other values of x\\ and x±: the agreement 
with the numerical results is very satisfactory as can be seen from Fig. 4. 
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We turn now to the double correlator 7° Eq. (7). We again use the non-abelian 
Stokes theorem and the cluster expansion, to represent 7 C as 



where the sum is over permutations of the order in which E\ and appear under 
the sign of the cumulant; the latter is denoted by double angular brackets and implies 
that vacuum insertions are subtracted from the averages of the field strengths. 
In our MC calculations partly reported in Table 2, the orientation of the plaquette P 
or P' was kept along the plane 14, and we should consider accordingly in (19) always 
Ei(x) or Ei(x') inside the cumulants. 

Now the symmetry requirements impose severe conditions on the nonzero values 
of ^ C nv,a'v'- Since both P and P' are chosen in the middle of the x 4 interval for 
the Wilson loop [— -|,-|], one can use the symmetry with respect to the change 
£4 — > —X4. In this way one can show that all odd-power averages of the type 
Ei(u)§Ei(y)§...Ei(P')& > should vanish, since they are odd with respect to 
£4 — > — £ 4 . 

This property holds if one inserts in this odd-power cumulant additionally several 
magnetic field operators B k ... 

Similarly, the correlator containing E>i(P) should depend on , e.g. 



and therefore should vanish whenever h k or hi are zero for both intervals, with = 
(u — x)^, — (u — £%. In this way one proves that (20) vanishes for i = 1 or 2 
identically, since h 3 = h' 3 = 0. For % — 3 the combination (20) should vanish when 
h 2 = h' 2 = 0. Using those criteria we keep in our Table 1 only results for 714,14, where 
the signal is largest. The contribution to 714,14 comes only from the quartic cumulants 
in (19), since as we discussed above, the triple correlator <^ E, QE^E^ 3> vanishes 
identically . 

One can conclude from the Table 1, that the quartic cumulant sharply decreases 
for large xj_, and the transverse shape of the string, using 7 C , is similar to that deduced 
from p c , the thickness of the string being of the order of 3 lattice spacings. 



l c ^'v'{x,x') = — J d yi dy 4 d Ul du 4 {^: E^yjQE^u^F^xjx (19) 
$F /i /„/(£ / )$ > +perm} 





< E^u^Biix^E^x')® >~ e ikl h k hi 



(20) 
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All the correlations 7 WM v , which should vanish on symmetry grounds, have been 
found to be zero within two standard deviations. In these higher statistics 

should be necessary for a clear answer. 

5 Conclusions 

One should stress three different aspects of the results reported above. First of all, 
we have presented the most detailed measurements of the connected correlators made 
so far. There is an agreement between the field distributions of this paper and those 
reported in [1], where only correlators of En were measured to be nonzero. Here 
with better statics all components of magnetic field B and electric field E are given 
for all possible configurations of one probing plaquette (for p c ) and two probing pla- 
quettes (for / y c ). Symmetry requirements impose severe restrictions on p c and 7 C 
(independently on cluster expansion) and predict zero results in most cases, expect 
for EniPu), ^2(^24) an d for j c (Eu, En). Results are compatible with the predictions, 
and hence statistical errors seem to be reliable. 

Secondly, our results serve to check the validity and usefulness of the cluster 
expansion for the signals like p c , 7 C . This expansion allows to express p c and 7 C 
in terms of a simpler (and more fundamental) bilocal correlator G 2 = (F<&F<&), 
which has been measured earlier [5] and thus yields a clear prediction for p c and 
7 C . Comparison with measured results made earlier in [6] and here in more detail 
in Fig. 3,4, shows a good agreement and supports the fundamentality of the bilocal 
correlator, which is known to define the nonperturbative dynamics of confinement [3]. 

In particular the asymptotics of the string profile at large x± is shown to be 
exponential, see eq. (18), just as the asymptotics of G 2 , measured in [5]. This in 
contrast to the behaviour inside the string, where the Gaussian-like flattening is 
observed before [1] and also in this paper. This effect is connected to the flattening of 
G 2 (x) at small x, necessary for its regularity at x = (Note that the latter property 
of G 2 was not taken into account in (15)) and also because of the smearing effect due 
to integration in (16), yielding polynomial factors in (18). 

Of special importance is the first estimate of quartic cumulant G 4 = ((F$) 4 ) 
through 7 C in (19). 

To obtain a more quantitative measure of G4, we represent the coordinate depen- 
dence of G 4 as an exponent similarly to G 2 (15) with the same p. In this case the 
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dimensionless ratio is 

gj(g) 7 C (0) 
(G 2 (0)) 2 P c (0))2 - 

This kind of estimate would result from the instanton-type vacuum with instantons 
of small size p, p < 0.3 /m and typical density of one instanton per /m 4 [7]. 

Finally, let us compare the transverse shape of the string (the string profile) mea- 
sured in Table 1 and analytically written in (18) with that of the dual Meissner 
effect. As is known from the theory of the type-II superconductors [8, Eq.(48.11)] 
the asymptotics of the magnetic field distribution off the vortex line is exponential, 
B{r) ~ exp(—r/5), where 5 is the so-called London's penetration length. One can 
see that the dual Meissner effect is capable of reproducing such fine details of the 
confinement picture, as the field distribution around the string. 

This conclusion agrees with a recent study made in a completely different approach 
in [9]. 
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Figure Caption 

• Figure 1: The operator p c . 

• Figure 2: The spatial structure of the parallel chromoelectric field E\\. 

• Figure 3: Different components of the field-strength tensor (Ej| and B», circles 
and diamonds respectively) vs. x\\ at fixed x± = 0. 

• Figure 4: E\\ vs. x±, for x\\ = 1,3,5,7 (triangles up, filled circles, diamonds 
and triangles down resp.). The dashed and solid lines are the results of the fits 
performed using resp. formula (18) and (3). 
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Figure 1: The operator p c . 
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Figure 2: Spatial distribution of E\ 
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Figure 3: E\\,B\\ vs. x\\, for x± = 0. 
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Figure 4: Eji vs. x±. 
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-.010907 ± 0.004644 
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5 


-.001269 ± 0.004222 


6 


0.000397 ± 0.003870 


Tab 


e 1: 7 C vs. x±, x» = 1. 
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